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Abstract. The existence of lower dimensional KAM tori is shown for a class of 
nearly integrable Hamiltonian systems where the second Melnikov's conditions are 
eliminated. As a consequence, it is proved that there exist many invariant tori and 
thus quasi-periodic solutions for nonlinear wave equations, Schrddinger equations 
and other equations of higher spatial dimension. 

1. Introduction and main results. 

Let us begin with the non-linear wave (NLW) equation 

utt - Uxx + V{x)u + h{x, u) = (1.1) 

subject to Dirichlet boundary conditions. The existence of solutions, periodic in 
time, for NLW equations has been studied by many authors. See [B, B-B, B-G,B- 
P, Br, L-S] and the references thcirin, for example. While finding quasi-periodic 
solutions, one will inevitably encounter so-called small divisor difficulty. The KAM 
(Kolmogorov-Arnold-Moser) theory is a very powerful tool to overcome the diffi- 
culty. This theory deals with the existence of invariant tori ( and thus quasi-periodic 
solutions ) for nearly integrable Hamiltonian systems. In order to obtain the quasi- 
periodic solutions of a partial differential equation, one may show the existence of 
the lower (finite) dimensional invariant tori for the infinitely dimensional Hamilton- 
ian system defined by the partial differential equation. Now consider a Hamiltonian 
of the form: 

i 

H ^ {uj,y) + ^njZjZj + Rix,y,z,z), l < oo, (1.2) 

with tangential frequency vector oj = (wi, ...,a;„) and normal frequency vector ft = 
(f2i, rij. When i? = 0, there is a trivial invariant torus x — cut, y ~ 0, z ~ z ^ 0. 
The KAM theory guarantees the persistence of the trivial invariant torus with 
a small deformation under a sufficiently small perturbation R, provided that the 
well-known Melnikov conditions are fulfilled: 

(fc, uj) - 7^ (the first Melnikov's) (1.3) 
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for all fc e Z" andl < i < l, and 

(fc, u) + ^1 - ^ (the second Melnikov's) (1.4) 

for all fc e Z" and 1 < 31,32 < hji 7^ .h- Sec [E,K1,P1,W] for the details. 
This KAM theorem can be applied to a wide array of Hamiltonian partial dif- 
ferential equations of 1- dimensional spatial variable, including (1.1). In that line, 
Kuksin[Kl,2] shows that there are many quasi-periodic solutions of (1.1), assuming 
that the potential V depends on an n-dimensional external parameter vector in some 
non-degenerate way. Wayne [W] obtains also the existence of the quasi-periodic so- 
lutions of (1.1), when the potential V is lying on the outside of the set of some 
"bad" potentials. In [W], the set of all potentials is given some Gaussian measure 
and then the set of "bad" potentials is of small measure. Bobenko & Kuksin[Bo-K] 
and P6schel[P2] get the existence of invariant tori and quasi-periodic solutions for 
a given potential ^(.t) = m G (0,oo). By the remark in [P2], the same result holds 
also true for the parameter values —1 < m < 0. When m G (—00, —1) \ Z, it is 
shown in [Yl] that there are many hyperbolic-elliptic invariant tori. More recently, 
the existence of invariant tori ( thus quasi-periodic solutions) of (1.1) are shown for 
any prescribed non-vanishing and smooth potential^ V{x) in [Y2] and for V{x) = 
in [Y3]. In [C-W], [Bol], [C-Y] and [Br-K-S], the equation (1.1) subject to periodic 
boundary conditions is investigated. 

For NLW equation (1.1) of spatial dimension 1, the multiplicity of normal fre- 
quency ilj is 1 in Dirichlet boundary condition and 2 in periodic boundary con- 
dition. Considering partial differential equations with spatial dimension> 1, a 
significant new problem arises due to the presence of clusters of normal frequen- 
cies of the Hamiltonian system defined by partial differential equations. In this 
case, the multiplicity of flj goes to 00 as \j\ 00; consequently, the second Mel- 
nikov's conditions is destroyed seriously, thus preventing the application of the 
KAM theorems mentioned above to Hamiltonian partial differential equations of 
higher spatial dimension. Bourgain[Bol-4] develops another profound approach, 
originally proposed by Craig- Wayne in [C-W], and successfully obtains the exis- 
tence of quasi-periodic solutions of the nonlinear Schrodinger (NLS) equations and 
NLW equations of higher dimension in space. Instead of KAM theory, the used ap- 
proach is based on a generalization of Lyapunov-Schmidt procedure and a technique 
by Frohlich and Spencer [F-S]. Bambusi calls this approach "C-W-B method". 

The advantage of the KAM approach is, on one hand, to possibly simplify the 
proof and, on the other hand, to allow the construction of local normal forms closed 
to the considered torus, which could be useful for the better understanding of the 
dynamics. For example, in generally, one can easy check the linear stability and the 
vanishing Lyapunov exponents. As a counterpart of the KAM theory, an advantage 
of the C-W-B method is more flexible than the KAM scheme to deal with resonant 
cases where the second Melnikov's conditions are violated seriously, although a 
(local) normal form can not be obtained by using the C-W-B method. Therefore, 
it is expected to find out a method by which the resonant cases can be easily dealt 
and by which the normal form can be obtained at the same time. In that direction, 
the first result is due to Bourgain[Bo5], to my knowledge. In order to introduce 
Bourgain's idea in [Bo5], let's recall the basic idea of lower dimensional KAM tori. 



'^This potential V contains no parameter. 
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See [E], [Kl] and [PI] for details. Consider a finitely dimensional Hamiltonian (1.2) 
with (. < oo. Decompose the perturbation i? in (1.2) into as follows 



R=R''{x) + {Ry{x),y) (1.5) 

+ {R^x),z) + {R'{x),z) (1.6) 

+ {R'%x)z, z) + {R''^x)z, z) + {R^^{x)z, z) (1.7) 

+ 0{\y\^ + \y\\z\ + \z\% (1.8) 



where [x, y, z, z) is in the some subset of T" x R" x x C, and where R^ , R^ : 
fn ^ C" and i?^, i?^ : T" ^ C and i?^^, R^^ are complex u x l matrices. In 
the classic KAM theory, one finds a series of symplectic transformations to kill the 
perturbing terms (1.5), (1.6) and (1.7) such that the transformed Hamiltonian is of 
the form 

H = {Uj,y)+Y,^jzjzj + 0{\y\^ + \y\\z\ + \z\'). (1.9) 

One sees easily that T" x {y}|j;=o x {z, z}z=z=o is an invariant torus of (1.9). Going 
back to (1.2), one can get an invariant torus of (1.2). The second Mclnikov's 
conditions (1.4) comes out while killing the perturbing term (1.7). Bourgain[Bo5] 
modifies essentially the idea mentioned above. He eliminates the perturbing terms 
(1.5) and (1.6) and puts the perturbing term (1.7) into the "integrable" part of the 
Hamiltonian H. Consequently, he obtains a Hamiltonian of the form 

H ={io,y) + Y,^jZjZj + {R''{x)z,z) + {R'%x)z,z) + {R''^x)z,z) 

j=i (1.10) 

+ 0{\yf + \y\\z\ + \zf). (1.11) 

This Hamiltonian H can be regarded as a counterpart of the local normal form 
of the classic KAM theory. It is clear that T" x {y = 0} x {z = z = 0} is still 
an invariant torus of H. Because of without killing (1.7), the second Melnikov's 
conditions (1.4) are not required. Therefore, using the Bourgain's idea, one can deal 
resonant cases where the second Melnikov's conditions (1.4) arc violated seriously; 
at the same time, one can get a normal form (1.10+11). However, the "integrable" 
term (1.10) is not really integrable, since it contains the angle- variable x in R^^{x)'s. 
Because of this fact, while eliminating (1.6), one will have to solve a homological 
equation with variable coefficients 

^/^{iJ, da:)F' + (f2° + R''{x) + ■ ■ ■)F' = R' (1.12) 

where is unknown function. By contrast, in the classical KAM theorem, the 
homological equation is of constant coefficients. To solve (1.12) one needs to inves- 
tigate the inverse of a "big" matrix A of the form: 

A= diag((A;,a;) + O°:fceZ",j = l,...,t) + (S^(fc-0 :A;,/eZ"). (1.13) 

As in the C-W-B method, Bourgain[Bo5] uses and develops the technique by 
Prohlich and Spencer [F-S] to investigate the inverse of A. Even if one gets the 
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inverse A^^ , one has to show that A^^ is off-diagonal decay, in order to control 
the exponent- weight norm of A~^. This technique depends heavily on the algebraic 
structure of {k,uj) + Qj. For nonlinear wave equation of spatial d-dimension, 

This makes the algebraic structure of (k,uj) + flj very intricate, if d > 1. 

In the present paper, we consider Hamiltonian (1.2) with i = oo. Following 
Bourgain's idea in [Bo5], we put (1.7) into the "integrable" part (1.10) and find a 
symplectic transformation to eliminate (1.6) and (1.5). This leads to solve the ho- 
mological equation (1.12) of variable coefficients. Thus, we also have to investigate 
the inverse of the "big" matrix A as in [Bo5] . However, we do not use the technique 
by Frohlich and Spencer [F-S]. We impose a symmetry condition on R as follows 

R{-x,y,z,z) = R{x,y,z,z) (1.14) 

where (x, y) is in some subset of T" x M" and z, z arc in some Hilbcrt space. With 
the symmetry condition (1.14) we can prove that the spectra of A are "twisted" 
with respect to u, thus that the spectra are different zero if digging out some w's 
of small Lcbesguc measure. This implies the existence of A~^. In addition, we 
do not need to know that A~-^ is off-diagonal decay, by choosing Sobolev space 
as our working space instead of an exponential weight space. Generally it is very 
hard to prove that A~^ is off-diagonal decay, which involves the intricate algebraic 
structure of {k,u!) + Qj. It should be noted that the symmetry condition (1.14) 
is fulfilled by many partial differential equations, such as nonlinear wave equation 
(1.1), nonlinear Schrodinger equation, nonlinear beam equation and others. 

The rest of the present paper is organized as follows. In §2., a KAM theorem 
to deal Hamiltonian (1.2) with t = oo is given out. This theorem does not require 
the second Afclnikov's conditions. In § 3., the KAM theorem is used to obtain 
invariant tori and quasi-periodic solution for nonlinear wave equation and nonlinear 
Schrodinger equation. The §4-7 are devoted to the proof of the KAM theorem. The 
fifth section is the essential part of the present paper. 

Acknowledgement. The present author is very grateful to Professor D. Bambusi 

for encouragement. 

2. A KAM theorem. 

2. 1 . Some notations. Denote by (^2 , 1 1 • 1 1 ) the usual space of the square summable 
sequences, and by (L^, || • ||) the space of the square integrable functions. By | • | 
the Euclidian norm. Let p > d/2. For a sequence z = (z^ e C : j G Z''), we define 
its norm as follows: 

INI^=Ebl'"l^il'- (2-1) 

Let P be the set of all sequences satisfying (2.1). It is easy to see that (P is a 
Hilbert space with an inner product corresponding to (2.1). Introduce the phase 
space: 

V := (C"/27rZ") x C" x F x (2.2) 
where n is a given positive integer. We endow V with a symplectic structure 

dx Ady + \/^dz Adz = dx Ady + ^ dzj A dzj , {x,y,z,z) e V. (2.3) 
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Let 

To" = (IR"/27rZ") X {y = 0} X {z = 0} X = 0} C -p. (2.4) 
Then 7q* is a torus in V. Introduce a complex neighborhoods of Tq in V: 

D{s, r) := {{x, y,u)&V: \lmx\ < s, \y\ < rM \z\\p < r, \ \z\ \p < r} (2.5) 

where r, s > arc constants. 

For r > we define the weighted phase norms 

,\W\p = \X\ + ^\Y\ + ^\\Z\\p + ^\\Z\\p (2.6) 

for W = {X, Y, Z, Z) G V . Let O c R" be compact and of positive Lebesgue 
measure. For a map W : D{s,r) x O ^ V, set 

f\W\pMB,r)xO ■■= sup f\W{x,0\p (2-7) 

and 

HW"lp,D(.,r)xo - max sup f\di,W{x,0\p, ^ = (Ci, ?«)• (2.8) 

Denote by C{£^, the set of all bounded linear operators from £p to f.^ and by 
III • I lip the operator norm. For any subset S C Z'^ with cardinality 5" < oo and a 
finitely dimensional vector u = {uj G C : j G S) and a matrix U = {Uij G C : i,j G 
S) of finite order, let 

u = {uj : j G Z'^), here Uj = 



Uj, j G S 
0, jGZ''\ S, 



and 



U = {Uij gZ'^), herei/y = 



[/y, i,j e S 
0, z or j e Z'' \ S. 

Define ||w||p := ||w||p and |||C/|||p := |||t/|||p. Similarly define ||u|| := ||u|| and 
infill •= infill- The following lemma states the relation between ||| • |||p and ||| • |||. 

Lemma 2.1. (i). Let Y = {Yij : \j\ < K) he K x K matrix with some K > Q. 
Then \\\Y\\\p < KP\\\Y\\\ and \\\Y\\\ < \\\Y\\\p. (it). Let Y = (Y,- : |j| < K) be a 
K -dimensional vector. Then \\Y\\p < KP\\Y\\ and \\Y\\ < \\Y\\p. 

Proof. The proofs are trivial. We omit it. 

In the whole of this paper, by C or c a universal constant, whose size may be 
different in different place. If / < Cg, we write this inequality as f < g when we 
dot not care the size of the constant C. Similarly, if / > Cg we write f > g. 

2.2. The statement of the KAM theorem. For two vectors b,c G C*^ or M*^, we 
write {b, c) — J2j=i ^j^j if < oo. If = oo, we write (&, c) = Y^^=i ^j^j- Consider 
an infinitely dimensional Hamiltonian in the parameter dependent normal form 

^0 = (c^°(e),y) + E {x,y,z,z) G V. (2.9) 
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The tangent frequencies ui^ = {cvf, • • • , w^) and the normal frequencies 0"'s (j e Z'^) 
depend on n parameters ^ £ Oq C M", where Oq a given compact set of positive 
Lebesgue measure. Let 

1^0(0 = diag(f}°(0 :jeZ''). (2.10) 

Let 

No = {LJ°{^),y) + {n%^)z,z), (2.11) 
where z = {zj : j G Z**). The Hamiltonian equation of motion of Nq are 

a; = w°(0, y = 0, i = \/^0°(0^, l=-\/^n°(0^- (2.12) 
Hence, for each ^ G Oq, there is an invariant n-dimensional torus 

To" = T" X {y = 0} X {z = z = 0} 

with frequencies ui^iO- The aim is to prove the persistence of the torus Tg*, for 
"most" (in the sense of Lebesgue measure) parameter vector ^ € Oq, under small 
perturbation R of the Hamiltonian iVo. To this end the following assumptions are 
required. 

Assumption A: (Multiplicity.) Assume that there are constants ci,C2 > such 
that for all ^ & Oq, 

(f^°)«<ci|jr^ (2.13) 
where we denote by (•)" the multiplicity of {■). 

Assumption B: (N on- degeneracy.) There are two absolute constant cs,C4 > 
such that 

sup \det d^Lu^iOl > C3, sup mu^Ol < Ci, j = 0,1. (2.15) 

Assumption C: (Analyticity of parameters.) Assume that both uj°{^) and fi"(^) 
are analytic in each entry S^i of the variable vector £^ = (^i, ^„) G Co and 
assume that both lo^(C) and r2°(C) are real for real argument ^. 
Assumption D: (Bounded conditions of Normal frequencies.) Assume that there 
exists constants C5,cq,C7 > and constant k > such that 

inf n° > csbT + C6, (2.16) 
€ee>o 

sup <C7«:l, Z = l,...,n (2.17) 

uniformly for all j. 

Assumption E: (Regularity.) Let so,rQ given. Assume the perturbation term 
R^{x, y, z, z; ^) which is defined on the domain D{so, ro) x Oq is analytic in the space 
coordinates and also analytic in each entry (I = 1, ...,n) of the parameter vector 

^ € Oq, and is real for real argument, as well as, for each ^ G Oq its Hamiltonian 
vector field X no := {Ry, —R'^, y/^dzR'^ , —^/--TdzR'^)^ defines a analytic map 

Xro : D(.sa,ra) cV 

where T = transpose. Also assume that Xj^o is analytic in each entry of ^ € Oq. 
Assumption F: (Symmetry.) For any {x,y, z,z,S,) € D{so,ro) x Oq, 

i?"(-x, y, z, z, = R"ix, 2/, z, z, 0- (2.18) 
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Theorem 2.1. Suppose H = Nq + satisfies assumptions A-F and smallness 
assumption: 

(2.19) 

Then, for < e <C 1, there is a subset Oc C Oq with 

Meas Oe > ( Meas Oo)(l - 0(e)), 

and there are a family of torus embedding ^ : x ^ V and a m,ap : Oe 
M" where and lo*{^) is analytic in each entry of the parameter vector 

^ = (^i,...,Cn) for other arguments fixed, such that for each ^ G O,, the m,ap $ 
restricted to T" x {^} is a, analytic embedding of a rational torus with frequencies 
'^*(0 /o*" the Hamiltonian H a,t ^. 

Each embedding is real analytic on x {^}, and 



1^ - ^o|p,TnxOe < ce, r-ol^ - ^olp.TnxO, < ce 



1/3 



w/iere $o i/ie trivial embedding T"' x Oq —^ T"' x {y = 0} x {z = z = 0}, 
c > 15 a constant depending on n. 

3. Application to nonlinear partial differential equations. 

3.1 Application to nonlinear wave equations. Consider nonlinear wave equation 

utt- Au + M„u + eu^ = 0, 0eT'',d>l (3.1) 
where u = u{t, 6) and A = is a real Fourier multiplier 

M„ cos(j, B) = a J cos(j, 9), M„ sin(j, 6) = sin(j, 6*), cr^ e M, j e Z''. (3.2) 

Pickaset e = {ei,...,e„} C Z''. Let Z^= Z'^\{ei, ...,e„}. Following Bourgain[Bo3] , 
we assume 

-e,=.„ (./ = l,...,n) 

= 0, j e Z^. 

Let (j G Z'') be the eigenvalues of the self-adjoint operator —A + subject 
to periodic b. c. ^ e T**, and let = (f>f{0) be the normalized eigenfunctions 
corresponding to . Then 

d 

>^f = bf + = T^ji + <^j, j = in, -Jd) e z^ (3.4) 
1=1 

and 

Note that {0^ : j e Z**} is complete orthogonal system in L'^iT'^). Let 

«(i,^)= E«.*W</'^^)- (3-6) 
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Inserting (3.6) into (3.1), then 

qf + Xfqf + e{u\ <^±)^2 = 0, j e Z'^. (3.7) 



Let g± = X Xfqf. Then 



^7 = ^^71 



3 -J 



(3.8) 



For simphfying notation, we do not distinguish + sign and — sign in (3.8). For 
example, we write as qj. However, we should keep in mind that qj runs over 
the set {qj',qj}- The system (3.8) is a Hamiltonian system with its Hamiltonian 
function 

H = H{q,q) = ^ ^ V^(«? + + ^G{q), (3.9) 



2 

3 



where 



and 



Let 



and 



^(9) = X] Gijkiqiqjquqi (3.10) 



Gijki = /, , ■ ■ : / (t)i(t>j(t^k<l>idx. (3.11) 



= = (i = l,...,n), (3.12) 



^° = V^i = VW, iezl. (3.13) 

Introduce a symplectic coordinate change 



Pei = + %)sina;j, qe^ = ^/2{l + yj)cosXj, j = 1, ...,n, 



To shorten notation, let 



yi = y/2{l +yi). 

Then (3.9) is changed into 



(3.14) 



H = H{x, y, z, z) = {to, y)+Y^ ^']zjZj + ei?° (3.15) 
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where 



■ ^ GijkiyiVjykyi cos cos Xj cos Xk COS Xi 

+ G ^ jkiy ty ]y k COS Xi COS Xj COS Xk — -1=- 

i,j,kee,leld 

+ 2^ Gijkmyj cosxicosxj — ^= K-^-^^) 

i,jee,kJelA ^ ^ 



Er^ ~ ^3 + + Zk Zi + Zi 
GijkWi COS Xi—^ -= ^ 



Zi + Zi Zj + Zk -\- Zk Zl -\- Zl 

ijkl " 



. ., , V2 s/2 s/2 

i,j,k,leM. V V V V 

It follows that in a suitable domain of (a;, y, z, z), 

R%-x,y,z,z)=R\x,y,z,z). (3.17) 

This implies the symmetry assumption F is fulfilled. By (3.5) and (3.11), it is not 
difficult to verify that 

G^Jkl = 0, imlcss i±j±k±l = 0, (3.18) 
for some combination of plus and minus signs. For w,z € the convolution z-kw 
is defined by [z ★ w)j = Y!,kei.'' '^i-k^k- 

Lemma 3.1. If p > d/2, then \\w-kz\\p < c||w||p||2;||p forw,z G with a constant 
c depending only on p. 

Proof. The proof is very elementary. See [p. 294-295, P2]. 

Let So = ro = 1, w° = {ui, a;„) and r2° = (r2° : j G Z^). And let parameter 
a = (fTi, cr„) runs over Oo ■= [1, 2]*^. Note that R^ is independent of a. By (3.18) 
and Lemma 3.1, we see that i?° is well defined on D{so, ro) and the Assumption E 
holds true; moreover, 

ro\XeR«\p.D(so,ro)xO < ro l-'^efl" lp,D(so ,ro) X O = 0- (3.19) 

This implies that (2.19) is fulfilled. It follows from (3.12) that the Assumption B 
(that is, (2.15) ) holds true. By (3.13), we see that the assumptions A, C, D are 
fulfilled. In particular, (2.16) holds true with k = 1. Using Theorem 2.1, we have 

the following theorem. 

Theorem 3.1. There is a subset Og c Co with Meas 0> {1 - Ce)Meas O such 
that for any a G Ojq, Eq.(3.1) with small e has a rotational invariant torus of 
frequency vector oj^ = a;°(cr). The motions on the torus are quasi-periodic with 

frequency . 

3.2 Application to nonlinear Schrddinger equations. Consider nonlinear Schrodinger 
equation 

^/^ut - Au + MaU + eu\u\^ = 0, 6* G T'', d > 1, 
where M„ is a Fourier multiplier 
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Theorem 3.2. Assume aj 's satisfy (3.3). There is a subset Oj^ C Oq with 

Measure O > (1 - Ce) Measure O 

such that for any a = (cti, cr„) e Oq, the nonlinear Schrddinger equation with 
small e has a rotational invariant torus of frequency vector w° = cJ^icr) where 

= le^r + j = l,-,n. 

The motions on the torus are quasi-periodic with frequency a;°. 

Proof. The proof is similar to that of Theorem 3.1. Note that in this case, (2.16) 
is fulfilled for k = 2. We omit the details. 

4. The linearized equation. 

4.1. Unperturbed linear system. Recall 0"'s satisfy the Assumptions A, C, D. 
Assume that there is a new tangent frequency vector^ a; satisfing the Assumptions 
B, C. Let O C Co be an open set. In this section, we pick two large constants K- 
and K. And let s = 1/K- and s' = A/K. (In the m-th KAM step, we will choose 
K_ = Km-i « 2(™-i)'/2 and K = Km- 2™V2.) Let 

^(s) = {a; G C"/27rZ" : < s}, s = l/K_. 

Let Xjv be a linear Hamiltonian system with Hamiltonian function: 

+ 1^2 mi{x;^>izj+ E m!i^^i)ziZi + l ^ ma^^ozizj) (4.i) 

ije'^'' i.jeZ'i i,jez<^ 

:={oj, y) + (0°^, z) + \{B^'z, z) + {B^'^z, z) + \{B!'z, z) 

where :B^j(a;;^)'s are analytic in a; G D{s) for any fixed ^ G O, and all of w(^), 
and Bij {x; ^) (for fixed x G D{s)) are analytic in each entry ^; (Z = 1, n) of ^ G O. 
Assume Bjj's satisfy the following conditions: 
(1.) Symmetry. 

m.^x, = B^Hx, 0, Bf/(-ar, = mf{x, (4.2) 

for any (x,^) G D{s) x O. 

(2.) Finiteness of Fourier modes. 

^?/(^;0= E (4-4) 

\k\<K 

^Here a; = a;(^) is not necessary to be up, due to the well-known phenomenon the shift of 
frequency in the KAM iteration. 
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where if 1 is a constant independent of x, ^, i and j. (In the m-th KAM iteration, 
wc will take K w 2'"'/2.) 
(3.) Boundedness. 

sup |||S(a;;Olli;<e«l, (4.5) 

D(s)xO 

where * = the blank and L. 

(4-) Non-degenerateness. Moreover, we assume 



ci > 



9e 



> C2 > 0, V ^ e O. (4.6) 



(This assumption enables us to regard u instead of ^ as a parameter vector.) 
2. split and estimate for small perturbation. We now consider a perturbation 

H = N + R 

where R = R{x,y,u;^) is a Hamiltonian defined on D{s,r) and depends on the 
parameter ^ g (D. Assume that the vector field : D(s,r) x O ^ V is real 
for real argument, analytic in {x, y, z, z) G -D(s, r) and in each entry of ^ € O. We 
assume that there arc quantities^ e = e{r,s,0) and e'-' = £^{r, s, O) which depend 
on r,s,0 such that 

r\XR\p,D{s,r)xO < r\Xjj\p jj(^g r)xO^^^' £ < < 1. (4.7) 

Let 

R= E E C(fc)e^('='^'y'"0''^^«^ (4.8) 

2|m| + |gi+g2|<2feeZ<i, 
meZ'*,gi,q2eZ~ 

with the Taylor- Fourier coefficients Rmqiq2{^) of ^ depending on ^ G O, and being 
analytic in each entry of 4- We see that is a partial Taylor-Fourier expansion 
of R. We will approximate R by R. Now we give some estimates of R. 

Lemma 4.1. 

r\XR\*p,D{s,r)xO < r-|-^Klp,r>(s,r)xO ^ S*' (4-9) 
nr\XR - ^^|*,i3(^,4^^)xO < V ■ r|^fl|*,£)(s,r)xO < V^* , (4-10) 

for any < r/ <C 1, where * = the blank or jC, for example, e* = e or e^. 
Proof. The proof is similar to that of formula (7) of [PI, 129]. 
With this lemma, we decompose R = R^ + R^ + R'^, where 

R'' = {R',z) + {R^,z), 

R^ = l{{R''z,z) + {R^'z,z} + {R''z,z}), 



^We will take e = e^''/^)"' and = s^/^ in the m-th KAM iteration step. 
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with i?^, Ry : D(s) X O ^ C"; i^^ : D{s) x O ^ IP; etc. Z)(s) x ^ 

£{iP,£P). For any vector or matrix 

we introduce the cut-off operator Tk as follows: 

rKY= J2 Y^g{k)e^^'''-'^y"'z''z^ (4.11) 

|;c|</f 

Lemma 4.2. Assume (s — s')ii' > | InTj]. Let Rk = R — ^kR- We have 

r\X(rKR)\*p,D{s',r)xO ^ r-|-'^fllp,_D(s,r) X O < ^* ' (4-12) 
r|^iiif|p,r>(s',r)xO < (4-13) 

where * = the blank or C. 

Proof. The proof of (4.12) is obvious. Let us give the proof of (4.13). Write 

Rk = R^k + iK^y) + {Rk, + {R]<, ^) + {Rk^, + {Rk'^, ^) + (RH^, z). 

Note that the terms Rj^, R^, and so on, are analytic in a; e D(s) for fixed ^ £ O. 
And observe that 

Rk= S|(fc)ev^('='-),-- - . 

\k\>K 

Then by Cauchy's formula, we have \R^{k)\ < e"*'*^' sup^j^^) l-R^^I, and so on. Thus, 

r\XRK\l,Dis',r)xO ^ r-|-'^Jj|p,£)(s,r)xO X] ^ V^* ■ 

\k\>K 

□ 

Now we can write 

R = rKR + Rk + {R- R)- (4.14) 
Let TZ = TkR- Hence we can write 

n =n^ + (7^^ y) + z) + (7^^ z) 

1 - - (4.15) 

+ -((7^^^^,^) + (7e^^^,^) + (7e^^^,^)). 

Noting (4.8) and (4.11) we have 

TV = Yj R*{k)e^^^'''\ (4.16) 

\k\<K 

where * = x,y, z, z, zz, zz, zz. 
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Lemma 4.3. Under the smallness assumption (4-V on R, the following estimates 
hold true: 

|a.7^^b(.)xo < Kr^s, |9.7^^|£(,)xo < Kr^^"^ (4-17) 

\ny\Dis).o < Ke, \ny\i^,)^a<Ke'' (4.18) 

mi\pM^)xO<Kre, \\n%^j,^^^^a < Kre"" , u&{z,z} (4.19) 

m''^\\\p,Dis)xO < Ke, m^''\\\inis)xO<Ke'', u,v e {z,z}. (4.20) 



Proof. Consider 7?.^^. Ohsevvethat R^^ = dzdzR\z=z=o- By the generalized Cauchy 
inequality (See Lemma A. 3 in [PI]), 

zz zz K 

lll^^^lllp,n(s)xO < ^lll^''''lllp,D(s)xO < —\\dzR\\p,D(s,r)y.O 

< r\XR\p^D{s,r)xoK < Ks. 

The remaining proof is simple. Wc omit the details. □ 

Lemma 4.4. Assume R satisfies the following symmetric condition: 

R{-x,y,z,z;0 = R{x,y,z,z;^), y {x,y, z,z;^) e D{s,r) x O. (4.22) 

Let Rfj , R^j , R^j be the elements of the matrices R^^,R^^,R^^, respectively. Then 
we have 

R%-x) = R'ix), Ry{-x) = Ry{x) (4.23) 

i?^(-x) = R\x), R%-x) = R^x) (4.24) 

i?-(x) = i?-(x), i?f/(-x) = R^fix) (4.25-1) 

R^fix) = R^fix), R^f{-x) = R^(x) (4.25-2) 

Proof. Noting that 

Rtj{x) = dz,dzjR{x, y, z, z; ^)|y=o,z=2=o, 

and dzidzj = dz^dz^ we get the first equation of (4.25-1). Applying dz^dz^ to both 
sides of (4.23), we get 

dzidz^ R{-x, y,z,z; ^)\y=Q^z=z=o 
=dzidzjR{x, y, z, z: C)|y=o,z=2=o 
=dzidz^R{x, y, z, z\ 0|y=o,z=2=o, 

that is, the second equation of (4.25-1) holds. The remaining proofs are similar to 
the previous one. □ 

Recall 7^ = VkR- Let 

^zz ^ Qzz _^ -ji^^ Q^^ = + 7^^^^ 522 = ^22 _^ 7^22 (4 26) 

uj = ui + Ry{Q). (4.27) 
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and 

N -. 

where 



(a;, y) + {n^z, z) + \{B^'z, z) + (B^^z, z) + z) 
-- B'^\ Then by (4.25), (4.2) and (4.3), 

i?-(x,e)=S|f(x,e), Sj7(-x,e)=i?f(a;,0 
(:r,e) =i3f (^,0, (-x,0 =Sf (^,0 



(4.28) 



(4.29) 



for any (a;,^) e £>(s) x £). Assume /Te < e and Ke'-' < e. By (4.5) and (4.20), we 

have 

<e (4.30) 

where * = the blank or C By (4.18) and Ke <C 1 as well as (4.6), we get that (4.6) 
still hold true for u: 



ci > 



det 



> C2 > 0, V ^ e O. 



In view of (4.27,28,14,15), 

H=N + Tl'' + {TV - Ry{0), y) + (7^^ z) + {U' , z) 
+ {R-TkR) + {R-R). 



(4.31) 



(4.32) 



4.3. DerivaMon of homological equations. The KAM theorem is proven by the 
usual Newton-type iteration procedure which involves an infinite sequence of sym- 
plectic coordinate changes. Each coordinate change is obtained as the time-1 map 
Xp\t^i of a Hamiltonian vector field Xp. Its generating Hamiltonian F solves the 
linearized equation 



{F, N} = TV' + {W - Rv{0), y) + (7^^ z) + (7^^ z) 



(4.33) 



where {•, •} is Poisson bracket with respect to the symplectic structure dx Ady + 

\J —\dz A dz. Without loss of generality, wc assume 7?.^(0) = since it docs not 
affect the dynamics. We are now in position to find a solution of the equation (4.33) 
and to give some estimates of the solution. To this end, we suppose that F is of 
the same form as the right hand of (4.33), that is, F = F^ -\-F^ , where 



F'' = F- + {F\y), 
F^ = {F\z) + {F-^,z), 



(4.34) 



with F'',Fv,F'', F^ depending on x, C, and F^(0) = 0, ^^(O) = 0. Let d^, = (w, 9^) 
where (•, •) is the usual inner product in M". Using (4.34) and (4.28) we can compute 
the Poisson Bracket {F,N}: 



{F,N}=-^{{Fy,d,)B''z,z) - {{Fy,d,)B''z,z) - ^((F^ a,)B-^, ^) 
+ d^F- + {d^Fy,y) + {d^F\ z) + z) 
+ ^/^((B^^F^ z) + {B'^^F^, z) + z)) 
rj^^^ff^z^^^ + (B'^^F^z) + (^°F^^)). 



(4.35) 
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Prom (4.33,35) we derive the homological equations: 

d^F- = TkR''{x,0, (4.36) 

d^py = TkR^x, - RKO), (4.37) 

-V^d^F'- +Vk{^'^F''- - B"-F'' + B'^'^F'') = -V^T kR"" (x, 0, (4.38-1) 

V^d^F' + TKin"F' - B"-F'- + B'-'^F'-) = ^/^T k R" [x , S) , (4.38-2) 

We should note that if F solves the equations (4.36-38), then F solves the following 
equation 

{F,N] =^{{Fy,d,B'')z,z) + {{Fy,d,B^'')z,z) + ^{{Fy,d,B^')z,z) 

+ V^{il-TK){B^'F'),z) +V^{{l-TK)iB''F'),z) (4-39) 
- V^((l - rK){B''F'),z) rK){B''F'),z) 

instead of (4.35). 

Lemma 4.5. If F in some sub-domain D{s', r') x C of D{s, r)xO is the unique of 
the homological equations (4-31-34), thenF satisfies the skew-symmetric conditions: 

F{-x, y, z, z; = -F{x, y, z, z;^, V {x, y, z, z; e D{s',r') x O'. (4.40) 

Proof By (4.29,30) and (4.38), we see that -F^(-x) solves (4.38-2). Thus, 

-F'{-x) = F%x) 

similarly, wc can show that F'^{—x) = —F^{x), Fy(—x) = —Py. Consequently 
F{-x, y, z, z; ^) = -F{x, y, z, z;^). □ 

4-4- Solutions of the homological equations. 

Proposition 1. ( Solution of (4.36).) There is a subset 0]_ C O with MeasO\ > 

(MeasO){l - R-^) such that for ^ G 0\, 

l(^''^(0)l>;^, for alio ^keZ^,\k\<K. (4.41) 

Then, on D{s') x 0]_, the equation (4-36) has a solution F^{x,^) which is analytic 
in X £ D{s') for ^ fixed and analytic in each ^j, {j = 1, ...,n) for other variables 
fixed, and which is real for real argument, such that 



I a 



-.nDisnxO^ < K'^r'e, < K'^r'e''. (4.42) 



Proof. The existence of the set 0+ is well-known in KAM theory. We omit the 
proof of the existence. Recall X-fi : D{s, r) cP ^ V is real analytic in {x, y, z, z) G 
D{s, r) and each entry of ^ € O. Expanding d^TZ^ into Fourier series 

0^11'== ^ a^(A:)e^('='^). (4.43) 
o#feeZ'',|fe|<A: 
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Since dxTV^ is analytic in x e D{s), we get that the Fourier coefficients 9a;7?.^(A;)'s 
decay exponentiaUy in k, that is, 

\d^{k)\ < \d,TZ-\nis)xoe-'^''^ < e-'^'^^Kr^s, (4.44) 
where we have used (4.17). Expanding S^-F^ into Fourier series: 

F"" = d7F^{k)e'^^^'''\ (4.45) 

07tfeez",|fe|<i<- 



Inserting (4.43,45) into (4.36), we get 



By (4.36,41) as well as Lemma A.l, we get that for x e D{s'), 

Applying d^. to both sides of (4.36) and using a method similar to the above, we 

can get the second inequality of (4.42). □ 

Proposition 2. ( Solution of (4.37).) On D{s') x 0\, the equation (4-37) has 
a solution F^^x,^) which is analytic in x G -D(s') for i; fixed and analytic in each 
= 1, ...,n) for other variables fixed, and which is real for real argument, such 

that 

Proof. The proof follows almost exactly the proof of Prop. 1. We omit it. □ 

5. Solution of (4.38). 

This section is essential part of the present paper. Let Bf^ e Z"^) be the 
elements of matrix B^^. Note that i?fj is a function of (x,^). Wc temporally omit 
the parameter ^ for simplifying notation. Regarding B^j as a function of x, we let 

Bl^{k) be the fc-Fourier coefficient of Bf--. Let 



Bt^^ = {BI^^{k-l):\k\,\l\<K). 

d: 

{i,jeZ'^), that is. 



Then Bf? is a matrix of order KK Let B be a block matrix whose elements are 



= {Bf? : i,j e Z'') = (Bffik - I) : \k\,\l\ < K;i,j e Z"^). 

Similarly we have B^^ and B^^ . Write F^ as a column vector of infinite dimension: 
F"" = {F^ e C : i e I/). Let Ff = {Ff{k) : \k\ < K) and F~^ = {fJ : j G U^). 
Similarly we have F^, W and W. Set 
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Introduce a tensor product space as follows 

P := F ® C^' ={z= {z, e C^^" : j € l") : (M;)' = ^ \z\'^ < oo}. 

j 

By abuse use of notation, we 1 1 • 1 1* as 1 1 • | |p. For any u,v G ff" , we can regard {u, v)'^ 
as a vector in ff' ® P'. Define 

||K^f||p = ^IHI, + ||^||,. 

Therefore, 

jez<i|fe|<if |fc|<is: 
where is fc-Fourier coefficient of F^{x). Similarly, 

11^11^= ^ u&{z,z}. 

\k\<K 



And 



In view of (4.29,30), we see that the linear operator B is self-adjoint in © i.e., 
B = B*. With these notations, we have the following lemma. 

Lemma 5.1. (i). The vector R is in ® with 

\\R\\p<Krs, \\R\\^<Krs'^ (5.1) 

(a). The linear operator B is self-adjoint in P ® with p = 0, i.e., B = B* ; and 

IPIIIp-III^IIIp <e (5-2) 

Proof. The conclusion (i) comes from (4.19). It follows form (4.30) that the con- 
clusion (ii) holds true. □ 

Let 

A± = diag (±(fc, w) +Qf'^:\k\<K,j&Z% K=(^ ^+ ^_ ) . 

Expand B^^.B^^, B^, F"^(x), F^{x), R''{x) and R^{x) in (4.38) into Fourier series, 
we get a lattice equation which reads 

{A + B)F = R. (5.3) 

Let 

M= (l + lOKm^lwl)^/". (5.4) 
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Write B = {Bij{kJ) : \k\, \l\ < K,iJ G Z^). Then we see that Bij{k,l) is one of 
B??{k - I), Blf{k - I) and Bl?{k - I). Let 



fill 


= {Bij{k,l) 


\k\,\l 


\<K,\i\ 


, \j\ < M) 




= (BijikJ) 


\k\,\i 


\<K,\i\ 


<M, \j\>M) 


B21 


= {Bij{k,l) 


\k\,\l 


1 < K, \i\ 


> M, \j\ < M) 


B22 


= {Bij{k,l) 


\k\,\l 


1 < K, \i\ 


> M, \j\ > M) 



Then 

And by (5.2), 
Let 

Then 

In view of (5.4), 



B 



\B, 



UMIP' I 



Bii B12 
B21 B22 



Biilllp <e, i,je{l,2}. 



(5.5) 



Ai = diag (±(fc, w) + : |fc| < K, \j\ < M) 
A2 = diag (±(/c, w) + 05 : |fc| < K, \j\ > M). 



A = 



Ai 
K2)' 



\±{k,uj) + ^°j\>K,\/\j\>M. 

Thus by (5.5) and (5.6), we see that there exists the inverse of A2 + B22 and 

00 

|||(A2 + B22)-^|||p< III E 111(^2 'S22)^lllpll|A2 'Hip <1- 
j=0 



Set 



B, 



lef 



Br 



El ^ 

-(yl2 + -B22)~^-B21 E2 

/ El — Bi2{A2 + B22)~^ 







Eo 



(5.6) 

(5.7) 

(5.8) 
(5.9) 



where Ei{E2, respectively) is a unit matrix of the same order as that of Ai (A2, 
respectively). It is easy to verify that 



\B 



Je/IIIP' lll-^HglllP 



< 1. 



Let 



Bii = Bii — Bi2{A2 + B22) ^B2i 
Then there exists the inverse oi A + B: 



{A + B)-'=B- 



i fiAi+Bii)-^ 



rig 







\ S-1 

{A2 + B22)-^l 



(5.10) 
(5.11) 

(5.12) 
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provided that there exists the inverse of Ai + Bu. And 

|||(^ + B)-i|||;< 111(^1 +^n)-^|||;, (5.13) 

where * = the blank or C. Now we are in position to investigate the inverse of 
Ai + Bii. First of all, we would like to point out that the matrix ISu is of finite 
order, and the order is bounded by 

K* := K^^'M'^ < ii:2n+d«-\ (5_;L4) 

Secondly, it follows from the sclf-adjoint-ness of B that the matrix Bu is also self- 
adjoint. Thirdly, by (5.5,7,11) we have 

IIIBiilll^ <e. (5.15) 

Fourthly, since each element of B is analytic in each entry of ^ G O, the matrix Bu 
is also analytic in each entry of ^ G O. Without loss of generality, we assume the 
first entry coi of lo is in the interval [1, 2]. Then the matrix Ai +Bii is non-singular 
if and only if ^Ai + lo^^Bu is non-singular, and 

|||(Ai +Bn)-'|||p < \\\{u^^A,+u>^^Bnr^\\\p. (5.16) 

Let 

S : ai = l/oJi, a2 = u)2/uJi,...,a„ = u>n/u)i. (5.17) 
Then it is easy to get 

d{ai, ...,a„) 



det 



d{u;i, w„) 



-(n+l) 



> 1. (5.18) 



Combining (5.18) and (4.31), we can regard a = (ai, ...,iJ„) as a parameter vector 
instead of ^, and we can show that 

Meas O < Meas S(0) < Meas O (5.19) 

as well as 

|P,Sn(?(a))||| <e, (5.20) 
d„a°j = 0(1), (5.21) 

where (2.17) is used in (5.21). Let 

n 

K = diag (±(/ei + km) + <yi^° ■ k = (fci, fc„) € Z", |A;| < K, |j| < M). 

(5.22) 

Take a as a parameter vector. Then 

wf^i +a;r^Bn(0 = K + <jSu{^{a)). (5.23) 
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Using Lemma A. 2, there are scalar functions /ii (cr) , ■ • • , /v a'* (o") and a matrix- valued 
function U{a) of order K* (See (4.14) for K*), which are analytic for each entry of 
real a and possess the following properties for every a G S(0): 

A, + aiBu{^{a)) =U{a)di&g{lii{a),--- ,HK'{^))U*{a), (5.24) 

and 

U{a){U{a)y = E. (5.25) 

It follows from (5.17) that" 

ll|t^(OIII = ll|f/1ll = l- (5.26) 

Arbitrarily take fi = fj,{a) € {/xi(ct), ...,^j.k* (c)}- Let 4> be the normalized eigenvec- 
tor corresponding jjL. Using Lemma A. 3 and noting (5.20) and (5.21), we get 

d,,li ={{d„M<^+aiBri{^{'j)))<l>A) 

=(diag(f20 : \j\ < M,\k\ < K)cj>,<P) + d,,{aiB{^{a)))cj>,<P) (5.27) 

= minO? + o(l). 

3 

Thus, 

9<TiM>c>0. (5.28) 

Let 

E{Oi) = {aeE{0) : < 1/{KK*)}, l = l,...,K*. (5.29) 
By (5.28), MeasS(0,) < 1/{KK*). Thus, 

K' 

Meas |JS(0;) < l/JsT. (5.31) 



Let 

Therefore, 

and for any a G S(0), 



E{0)=E{0)\[jEiO,). (5.32) 
1=1 



Meas E{0) > Meas S(0)(1 - 0(-^)) (5.33) 



|w(a)| > 1/KK*. (5.34) 

Moreover, 

IIRA,, + aiBn)-^||| < max{|Mr'|} < KK* < (5.35) 

where C = n + dK~^ + 1. Returning to the parameter ^, by (5.19), there is a subset 
OlcO with 

Meas Ol > (Meas 0)(1 - CiK'^), (5.36) 



"Here 1 1 1 • 1 1 1 is the ^2 norm of matrix. 
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and for any ^ e 

||i(Ai+Bn)-i||| <C2i^^. (5.37) 

Observe that /i; = is continuous in ^. By means of adjusting the constants 

Ci,C2 in (5.36,37), we can extend the set to an open set such that (5.36,37) 
still hold true for any ^ in the open set. Still denote by the open set. Replacing 
the i?2 norm ||| • ||| by ||| • |||p in (5.37), in view of the order of the matrix Ai + Bn 
being bounded by K* , we get 

|||(Ai +Sn)-i|||p < K^{K*)P < (5.38) 

where C is a constant depending n,d,K,p. By (5.13), 

\\\{A + B)-^\\\p<K^,^€Ol. (5.39) 

Note \\\d^{Ai + B)\\\p < K. We have 

\\\d^{K + B)-^\\l<K^+^, ^&Ol. (5.40) 

By (5.1) and (5.39,40), 

\\F\\p<K'^re, \\d^F\\j, < K^^ re^ , (5.41) 
where (7 is a constant depending p, d, n, k. Recall 

mi= E + E WF'mi- (5-42) 

\k\<K \k\<K 

Recall 

s' = -i. (5.43) 

Then for u G {z, z}, 



i;)||2 
\p 



That is. 



sup \\F"{x,0\\l= sup II E 

<K J2 ii^(fc)iiy'=i"' 

\k\<K (5.44) 

\k\<K 

<e'K\\F\\l<{K^rsf 



sup ||F"(a;,^)||p < K^rs, u € {z,z}. (5.45) 

D{s')xOl 



Similarly, 

Consequently, we have the following proposition. 



sup \\d^F^{x,0\\p<K''re, ue{z,z}. (5.46) 

D{s')xOl 
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Proposition 3. There is a subset 0^(10 with Meas 0% = {Meas 0){1-0{1/K)) 
and there is functions F^,F^ : D{s') x which are analytic in x G -D(s') 

and also analytic in each entry of £^ £ 0\; moreover the functions F^,F^ solve 
(4-38) and satisfy the estimates (5.45,46). 

6. Symplectic change of variables. 

In this section, our procedure is standard and almost tlie same as that of Section 
3 in [PI, p. 128-132]. Here we give out the outline of the procedure. See [PI] for the 
details. In this section, we denote by C a universal constant depending only on 
n,d,p,K. The constant C might be different in different places. 

6.1. Coordinate transformation. By Propositions 1-3, wo get a Hamiltonian F 
on D{s', r) where 

F = + {Fy,y) + {F',z) + (F^ z) 

and give estimates of F^, F^, F^ and F^. Let Xphe the vector field corresponding 
to the Hamiltonian F, that is, 

Xf = i-dyF, d^F, V^d^F, 

here and dz are the usual £^-gradients. Recall s = A/K- and s' = A/K. And 
let s" = 2,IK, s'" = 2/K, s"" = 1/K. Let 

It follows from Prop. 1,2 and 3 that for {x,y,z,z;^) G D{s',r)x G 0+, 
AXfI =\dyF\ + ^\d.F\ + + l\\d,F\\, 

=\Fy\ + ^\d.F'^\ + ln + l\\F% 

<\Fy\ + ^\F^\ + ^\\F'\\p + ^\\F'\\p 

where C is a large constant depending only on n, k, d and p. That is, 

r\XF\p,D{s',T)xO+ < K'-'e. (6.1) 

Similarly, we have 

r\XF\p^D(s',r)xO+ < K'^ ■ (6.2) 

where we have used the assumption e < e'-'. As in [PI, p. 129], we introduce the 
operator norm 

r\\L\l=snp^^^. (6.3) 

Using (6.1), (6.2) and the generalized Cauchy's inequality (See Lemma A. 3 of 
[Pl,p.l47]) and the observation that every point in D{s",r/2) has at least r\ ■ \p 
distance 1/9K to the boundary of D{s' , r), wc get 

sup r||Z)Xi?||p < K r\XF\p,D(s' ,r)xO+ 

< K^e. (6.4) 

D(s",r/2)xO+ 
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sup A\DXf\\^<K ,\XF\^^^j,^^,^,^^a^<K^e'', (6.5) 

_D(s",r/2)xO+ 

where DXp is the differential of Xp. Assume that K'-^ e and K'-'e^ are small 
enough. (These assumptions will be fulfilled in the following KAM iterations. In 
fact, in m-th KAM step, K « C2"'/2 ^nd e = e^^/s)". it follows that K^e < 1.) 
Arbitrarily fix ^ G 0+. By (6.1), the flow Xp of the vector field Xp exists on 
D{s"',r/4:) for t G [—1,1] and takes the domain into D(s",r/2), and by Lemma 
A.4 of [PI, p.l47], we have 

r\Xp - id\p^D{s"',r/4)xO+ < r\Xp\p^D{s' ,r)xO+ < K*-^ ^ (6.6) 

^^'^ t C 

r\Xp — «rf|j,_£)(s"',r-/4)xO+ 

< e^V{r\\DXp\\p^D{s",r/2)xO+) ■ r-|-^F|p,D(s',r)xO+ (6-''') 

for t e [—1, 1]. Furthermore, by the generalized Cauchy's inequality, 

r\\DX*p - I\\p,D(s"",r/8)xO+ < K^S, t G [0, 1] (6.8) 

and 

r\\DX*p - I\\inis"",r/s)xO+ < K"^^", t e [0, 1] (6.9) 

6.2. The new error term. Subjecting H = N + R to the symplectic transfor- 
mation $ = X^|t=i we get the new Hamiltonian scale := H o ^ = H o Xp on 
D{s"" , rjr) where < 77 -C 1. By Taylor's formula 

H+ = {N + R) o Xj. = {N + R + {R - R)) o X^, 
= {N + -R + Rk + {R-R))oX^p 

= N-{F,N}+ {t{F,N},F}oX*pdt (6-10) 
Jo 

+ -R+ [ {']l,F}oX'p dt + {RK + {R-R))oXl,. 
Jo 

Recall that (4.39) holds true when F solves (4.36-38). Thus, 

H+=N+ + R+ (6.11) 

where 

N+ = N+ i((F^ a,S-)z, z) + {{Fy, d,B^')z, z) + i((F^ 9,5^"^")^, z) (6.12) 

R+ = R\ + R\ + R% (6.13) 

where 

R\ =^{{l - Tk){B-^-^F^Iz) + ^/^((l - Tk){B^-^F^)z) 

_ _ , _ _ (0.14) 

1 ((1 - Tk){B''F~'), z) Tk){B''F'), z) 
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Rl = RKoX}r + {R-R)o Xjr (5.15) 



Rl= f {{1 + t)TZ, F} o X'p dt. 
Jo 

By (4.26,27,28), 



(6.16) 



with 



N+ = ico+, y) + {Q'^z, z) + \{Bl^z, z) + {Bfz, z) + \{Bfz, z) (6.17) 

= + i?^^ + {Fy,d^iBf^ + i?^^)), (6.18) 

^zz ^ QZ-z ^ ^zz ^ ^py^ Q^^^zz ^ ^zz-^-^^ 



and 



uj+=iO = w + Fy{Q). (6.19) 
Hence, the new perturbing vector field is 

= ^k^ + (^F)*(^fl -Xr + Rk)+ j (X^)*[X(i+t)^,Xf] dt, 

+ Jo 

where {Xyy is the pull-back of X^ and [•,•] is the commutator of vector fields. 
We arc now in position to estimate the new perturbing vector field Xp^^ . Let Y : 
D{s' ,r) C P — > "P be a vector field on D{s' , r), depending on the parameter ^ G 0+. 
Let U = D{s"",r]r) x 0+ and V = D{s'",2r)r) x 0+ and W = D{s",4r)r) x 0+. 
By (6.6) and the "proof of estimate (12)" of [PI, p.131-132]^, we have that 

r,r\{X'p)*Y\p,u < vr\y\p,V (6.20) 

and 

K 

vAiXpyYl^^ < nr\Y\py + nr\Y\p,W ■ r,r\XF\py. (6.21) 

We assume that 

eK'^/tj^ < 1. (6.22) 

These assumptions will be fulfilled in the KAM iterative lemma later. (In fact, in 
m-th KAM iterative step, we will let K'^ w 2<^™^£: = e'^'" ,r] = e^l'^. This implies 
(6.22) is fulfilled.) By (4.10) and (6.20,21), 

^r\{X\Y{X^-Xn%,v < rir\Xi,-Xn\py <r]e (6.23) 

and 

^r\{X\r{X^ - Xn)\iv < ne" + -^e" < r^e" . (6.24) 

Recall that (4.9) holds still true after replacing R by 7?.. By (4.2) and (5.4,5) and 
using the generalized Cauchy estimate, following [PI, p. 130-131] we get 

r\\X-'R,(t),XF^p^U < K r\Xn\p,V ' rl-'^Flp.y 

< K r\X^\p^w ■ r\XF\p,W (6.25) 



^Let a = in (12) of [PI]. 
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and 

r\[Xn(t),XF]\p.u 

< K r\Xf^\p r\XF\p,W + K rl-'^fllp.W r|-''^F|p,VF (6.26) 

< K^e^e + K^ee^ < rje^ 

Finally, we have 

r,r\Y\p,u < r,-^ r\Y\p,U, r,r\Y\iu < V'^ r\Y\§,u, (6-27) 

for any vector field Y. Note that for any function / : D{s) x O — > or M" which 
is analytic in a; e D{s) and in each entry of ^ e O, 

\\{l-rK)f\\Dis')xO< E 

\k\>K 

<||/lb(.)xO E e-(-^')l'=l (6.28) 



\k\>K 



< 



D(s)xO ■ ' 



.2 



where we have used 

(s-s')^~2|ln£|. (6.29) 
Applying (6.27,28) to (6.14), we can easily get 

rir\XfiiJp,u < ??"^e^ < V£:, vAXfiiJp^u < V^^- (6-30) 
Collecting all terms above, we then arrive at the estimates 

r+\XfiJp,D{s+,r+}xO+ < V^, r+ |p,_D(s+ ,r+) X 0+ < '7^'^ • (6.31) 

where 

s+ = s"", r+ = rjr. (6.32) 

It follows from (4.22) and (4.40) that the perturbing Hamiltonian R+ satisfies the 
same symmetric condition as R, that is: 

R+{-x,y,z,z;^) = R+{x,y,z,z;^), \/ {x,y, z,z;£,) e D{s+,r+) x 0+. (6.33) 

Finally we can easily check that 7V+ satisfies the same conditions as B_, that is, the 
conditions (4.2-6) hold true by replacing _B by 5+ and replacing w by w+. Here we 
omit the details. 

7. Iterative lemma and proof of the theorem. 

7.1. Iterative constants. As usual, the KAM theorem is proven by the Newton- type 
iteration procedure which involves an infinite sequence of coordinate changes. In 
order to make our iteration procedure run, we need the following iterative constants: 

1. eo = e, Q=e^(4/3)',/ = l,2,...; 

2. 7y; = e,'/y = 0,1,2,..; 

3. Ki = 2S^=iJ|lne|, I = 1,2,..., (thus, Ki = 2'('+i)/2| ine|); 

4. So > is given, si = l/if;-i, 1 = 1,2, (thus, (s; - si+i)Ki > 2| lne;|.); 
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5. To > is given, r; = r/,ro, ^ = 1, 2, 

6. £1(50 = {xe C"/(27r)" : |3a;| < s/} 

7. £1(5,, n) = {(x,y,z,z) e V : \Qx\ < .3i,\y\ < rf,\\z\\i < ri,\\z\\p < n}. 

7.2. Iterative Lemma. Consider a family of Hamiltonian functions Hi (0 < I < m): 



+ \{Bf\x, i)z, z) + Ri{x, y, z, z, ^ 



(7.1) 



where B^^ {x; S^) 's are analytic in x G D{si) for any fixed G Oi, and all of 
u}i{£,),n^{^) and Bi^- {x: S,) 's (for fixed x G D{si)) are analytic in each entry 
(I = 1, ...,n) of € Oi, and Ri{x, y,z,z,S,) is analytic in {x, y,z,z) G D{si,ri) and 
analytic in each entry of ^ G Oi. Assume Bf^ = {B^j : i,j G Z'*) '5 satisfy the 
following conditions: 

(1.1) . Symmetry. For any (x,^) G D{si) x Oi, the operators Bi^.,B^^ and B^z are 
self-adjoint from li to (.2- 

(1.2) . Finiteness of Fourier modes. 

\k\<Ki 



(1.3). Boundedness. 



sup \\\Brix;m\;<e, ... 

D{s,)xO, 



where * = the blank and L. Moreover, we assume the parameter sets Oi's satisfy 
(U). 

with 

Meas Oi > ( Meas Oo){l - K^^)] 
{1.5). The map ^ 1— > is analytic in each entry of ^ G Oi, and 



inf 

Oi 



det 



> ci > 0, sup \diu)i\ < C2, j = 0, 1. 

Oi 



{1.6). The perturbation Ri{x,y,z,z;^) is analytic in the space coordinate domain 
D{si,ri) and also analytic in each entry (k = l,...,n) of the parameter vector 

^ G Oi, and is real for real argument; moreover, its Hamiltonian vector field X := 

{dyRi,—dxRi,\/—ldzRi,—\/^dzRi)^ defines on D{si,ri) a analytic map 

Xf,^:D{si,ri)cV^V. 

{1.7). In addition, the vector field X^^^ is analytic in the domain D{si,ri) with small 
norms 
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Then there is is an absolute positive constant e* enough sm,all such that, if < 
€o < e* , there is a set Om+i C Om, and a change of variables ^m+i '■ T^m+i •= 
D{sm+i,rm+i) X Ora+1 ^ -D(sm, r^) being real'', analytic in 

{x,y,z,z) € D{Sm+l,rm+l) 

and each entry ^ G Om+i, as well as following estimates holds true: 

and 

Furthermore, the new Hamiltonian Hm+i ■= o ^m+i of the form 

satisfies all the above conditions {l.l — 7) with I being replaced by m + 1. 

7.3. Proof of The Iterative Lemma. 

As stated as in the iterative lemma, we have got a family of Hamiltonian functions 
Hi's {I = 0, l,...,m) which satisfy the conditions {l.l — 7). We now consider the 
Hamiltonian Hm- Let N_ = N^, R = Rm, ^+ = -^m+i and R+ = Rm+i- Let s 
- s"" = Sm+i, V = Vm, r = rm = VrnrQ, s = e,„, = eK^ and Om+i = 0\U0\. 
Clearly, e < e^. By means of the conclusion in Section 4, we got that there is a 
Hamiltonian F = Fm defined on'' D{sm+i,'fm+i) x Cm+i and a symplectic change 
of variables <I>m+i = Xp^\t=i. Note r;„e„j = em+i. It is easy to verify that the 
conditions (m + 1.1 — 7) are fulfilled. We omit the details. 

7. Proof of the theorem 2. 1 

The proof is similar to that of [PI]. Here we give an outline. By Assumptions 
A,B, C, D, E and the smallness assumption in Theorem 2.1, the conditions (LI — 7) 
in the iterative lemma in Section 6.2 are fulfilled with I = 0. Hence the iterative 
lemma applies to H. Inductively, we get what as follows: 

(i) . Domains: for m = 0, 1, 2, 

(ii) . Coordinate changes: 

^f™ = $1 o . . . o : V^+x D{so, ro), ; 

(iii) . Hamiltonian functions Hm {in = 0,1,...) satisfy the conditions (LI, 2, 3) with I 

replaced by m; 

®The word "real" means $m+i(z,^) = for any (z,{) G T>m+i- 

^Note Vm+i ■■= D{ ) X Om+1 C D{s,n+l,rm+l) X Om+l- 
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Let Oao = n^^gOm, 2?oo = ^'Dm- By the same argument as in [PI, pp.134], we 
conclude that ^"^,D^"^,Hm,^Hm converges uniformly on the domain I>oo, and 
Xjj^ ° = -0*°° • X^, where 

:= lim = (MO^V) + {n°{Oz,z)+ 

m^oo 

+ ^{B^^{x,Oz,z) + ^{B^J{x,Oz,z) + ^{B'J{x,Oz,z) 

here = lim^^oo B^,..., and X,^^ is the constant vector field on the torus T". 
Thus, X {0} X {0} is an embedding torus with rotational frequencies e 
Lu^{Ooc) of the Hamiltonian i^oo- Returning the original Hamiltonian H, it has 
an embedding torus $°°(T" x {0} x {0}) with frequencies This proves the 

Theorem. □ 

9. Appendix A. Some Technical lemmas. 

Lemma A.l. For fi > 0,1/ > 0,the following inequality holds true: 

^e-^l'=l'^|fcr<(^r-^(l + e)''. 
Proof. This Lemma can be found in [B-M-S] . 

Lemma A. 2. Consider an n x n complex matrix function Y(^) which depends on 
the real parameter ^ G M. Let Y{^) be a matrix function satisfying conditions: 

(i) Y{£) is self-adjoint for every ^ G M; i.e., Y{^) = {Y{^))*, where star denotes the 
conjugate transpose matrix; 

(a) l^(^) is an analytic function of the real variable ^. 
Then there exist scalar functions ^''^'^ ^ matrix-valued function 

U{^), which are analytic for real ^ and possess the following properties for every 
^ G M; 

Y{0 = U{Odiag{,,^{0,--- ,l^n{0)U*{0, U{0{U{Or = E. 

Proof See [pp.394-396, G-L-R]. □ 

It is worth to point out that this lemma does not hold true for ^ e M*^ with 
k>l. See [Ka]. 

Lemma A. 3. Assume Y = F(^) satisfies the conditions in Lemma A. 3. Let fx = 
/Lt(^) be any eigenvalue of Y and cf) be the normalized eigenfunction corresponding 
to 11. Then 

d^ii = {{d^Y)ct>, (!>)). 

Proof. The proof can be found in [Ka,p.l25]. 
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